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Abstract 

A search for the decay D*^7r^ is presented using a data sample corresponding 

to an integrated luminosity of 1.0 fb~^ ofpp collisions collected by LHCb. This decay 
is expected to be mediated by a H^-exchange diagram, with little contribution from 
rescattering processes, and therefore a measurement of the branching fraction will 
help to understand the mechanism behind related decays such as — )• tt+tt" and 
Bg DD. Systematic uncertainties are minimised by using B^ D*^tt^ as a 
normalisation channel. We find no evidence for a signal, and set an upper limit 
on the branching fraction of B{B^, D*^tt^) < 6.1 (7.8) x 10"^ at 90% (95%) 
confidence level. 
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Decays of mesons to final states such as D^D^ , D^D^ [T] and tt+tt" [2] liave been 
recently observed by LHCb. Such decays can proceed, at short distance, by two types of 
amplitudes, referred to as weak exchange and penguin annihilation. Example diagrams 
are shown in Fig. [TJ^a) and (b). There is also a potential long distance contribution 
from rescattering. For example, the D^D~ final state can be obtained from a 6 — )■ ccs 
decay to DfD~ followed by the ss pair rearranging to dd. Understanding rescattering 
effects in hadronic B meson decays is important in order to interpret various CP- violating 
observables. 




(a) (b) (c) 

Figure 1: Decay diagrams for (a) 5° via weak exchange, (b) 

via penguin annihilation, (c) — > D^*^~7t'^ via weak exchange. 



A measurement of the branching fraction of the decay B^ D*^7i^ can be used to 
disentangle the contributions from different decay diagrams and from rescattering [3lll]. 
This decay has only weak exchange contributions, as shown in Fig. [T](c). (The suppressed 
diagram for B^ — D*~^7r~ is not shown.) Moreover, rescattering contributions to the 
5^ — >■ D^*^^Tr'^ decay are expected to be small [S]. Therefore, if the observed branching 
fraction for the decay B^ ti'^ti" is explained by rescattering, a low value of i5(-B° 
D*^7i^) = (1.2 ± 0.2) X 10^^ is predicted [5]. However, if short-distance amplitudes are 
the dominant effect in B^ n^n~ and related decays, B{B^ D*~7i^) could be much 
larger. The measured B^ — > DD [1] and B^ — > |6] rates are at the upper end of the 
expected range in the rescattering-based model, but further measurements are needed to 
establish whether long-distance processes are dominant in these hadronic B decays. 

In this paper, the result of a search for the decay B^ -^^ D*^7i^ is presented. No 
previous measurements of this decay have been made. The inclusion of charge conjugated 
processes is implied throughout the paper. Since the flavour of the B^ meson at production 
is not tagged, the D*~7r~^ and D*^tt~ final states are combined. The analysis is based on 
a data sample corresponding to an integrated luminosity of 1.0 fb~^ of LHC pp collision 
data, at a centre-of-mass energy of 7TeV, collected with the LHCb detector during 2011. 
In high energy pp collisions all h hadron species are produced, so the B^ — t- D*^^^ decay, 
with branching fraction B{B^ D*-7r+) = (2.76 ±0.13) x lO'^ [71|8], is both a potentially 
serious background channel as well as the ideal normalisation mode for the measurement 
of the B^ branching fraction. 

The LHCb detector [S] is a single-arm forward spectrometer covering the 
pseudorapidity range 2 < 77 < 5, designed for the study of particles containing 6 or c 
quarks. The detector includes a high precision tracking system consisting of a silicon-strip 
vertex detector surrounding the pp interaction region, a large-area silicon-strip detector 
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located upstream of a dipole magnet with a bending power of about 4 Tm, and three sta- 
tions of sihcon-strip detectors and straw drift tubes placed downstream. The combined 
tracking system has momentum resolution Ap/p that varies from 0.4 % at 5 GeV/c to 0.6 % 
at 100 GeV/ c, and impact parameter (IP) resolution of 20 |J.m for tracks with high trans- 
verse momentum (px)- Charged hadrons are identified using two ring-imaging Cherenkov 
detectors. Photon, electron and hadron candidates are identified by a calorimeter system 
consisting of scintillating-pad and preshower detectors, an electromagnetic calorimeter 
and a hadronic calorimeter. Muons are identified by a system composed of alternating 
layers of iron and multiwire proportional chambers. 

The trigger [10] consists of a hardware stage, based on information from the calorimeter 
and muon systems, followed by a software stage which applies a full event reconstruction. 
In this analysis, signal candidates are accepted if one of the final state particles created 
a cluster in the calorimeter with sufficient transverse energy to fire the hardware trigger. 
Events that are triggered at the hardware level by another particle in the pp hbX 
event are also retained. The software trigger requires characteristic signatures of 6-hadron 
decays: at least one track, with > 1-7 GeV/c and Xip with respect to any primary 
interaction vertex (PV) greater than 16, that subsequently forms a two-, three- or four- 
track secondary vertex with a high sum of the of the tracks and significant displacement 
from the PV. The Xip is the difference between the of the PV reconstruction with and 
without the considered track. In the offline analysis, the software trigger decision is 
required to be due to the candidate signal decay. 

Candidates that are consistent with the decay chain B^^-^ Z}*^7r^, D*^ — ^ D^n^, 

A'+TT^ are selected. The and D*^ candidate invariant masses are required to sat- 
isfy 1814 < TUK+n- < 1914MeV/c2 and 2008.78 < m^o^- < 2011.78 MeV/c^, respectively, 
where a mass constraint is applied in the evaluation of m^o^-. The bachelor pion, 
from the B^^^ decay, is required to be consistent with the pion mass hypothesis, based 
on particle identification (PID) information from the RICH detectors [H]. All other se- 
lection criteria were tuned on the B^ — >■ D*^7r^ control channel in a similar manner to 
that used in another recent LHCb publication [T^]. The large yield in the normalisation 
sample allows the selection to be based on data, though the efficiencies are determined 
using Monte Carlo (MC) simulated events in which pp collisions are generated using 
Pythia 6.4 [13] with a specific LHCb configuration [H]. Decays of hadronic particles are 
described by EvtGen [15]. The interaction of the generated particles with the detector 
and its response are implemented using the Geant4 toolkit [IS] as described in Ref. [T7] . 

The selection requirements include criteria on the quality of the tracks forming the 
signal candidate, their p, px and inconsistency with the hypothesis of originating from 
the PV (xip)- Requirements are also placed on the corresponding variables for candidate 
composite particles B'^g)) together with restrictions of the decay fit (Xvertex); the fiight 
distance (xiight)' the cosine of the angle between the momentum vector and the line 
joining the PV to the 5°^-) vertex (cosOdh) [18] . 

Further discrimination between signal and background categories is achieved by calcu- 
lating weights for the remaining B'^ candidates [TP]. The weights are based on a simplified 
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fit to the B candidate invariant mass distribution, where the 5° region is neither exam- 
ined nor included in the fit. The weights are used to train a neural network [20] in order 
to maximise the separation between categories. To retain sufficient background events for 
the network training, the requirement on m^o^- is not applied. A total of fifteen vari- 
ables are used as input to the network. They include the Xip of the four candidate tracks, 
the Xip, Xvcrtex' Xflight cos6'dir of the and 5°^-) candidates, and the -8°^,^ candidate 
Pt. The pt asymmetry and track multiplicity in a cone with half- angle of 1.5 units in 
the plane of pseudorapidity and azimuthal angle (measured in radians) [21] around the 
B^^^ candidate flight direction are also used. The input quantities to the neural network 
depend only weakly on the kinematics of the B^^^ decay. A requirement on the network 
output is imposed that reduces the combinatorial background by an order of magnitude 
while retaining about 75 % of the signal. Potential biases from this data-driven method 
are investigated by training the neural network with different fractions of the data sample. 
The same results are obtained using a neural network trained on 30, 40, 50, 60 and 70 % 
of the total data sample. 

After all selection requirements are applied, approximately 50 000 candidates are se- 
lected in the invariant mass range 5150 < ■m£,,-j^+ < 5600MeV/c^. About 1% of events 
with at least one candidate also contain a second candidate. Such multiple candidates 
are retained and treated the same as other candidates. 

In addition to combinatorial background, candidates may be formed from misiden- 
tified or partially reconstructed 5°^^ decays. Contributions from partially reconstructed 
decays are reduced by requiring the invariant mass of the B^^^ candidate to be above 
5150MeV/c^. The contribution from B^^^ decays to identical final states but without in- 
termediate charmed mesons is negligible due to the requirement on the D*~ candidate 
invariant mass. A small but significant number of background events are expected from 
B^ — )• D*~K^ decays with the misidentified as a pion. The branching fractions of 
5° — )■ D*~K^ and D*~p are expected to be small due to CKM suppression, so that 

these potential backgrounds are negligible. 

Since the B^ decay mode is several orders of magnitude more abundant than the 
B^ decay, it is critical to understand precisely the shape of the B^ signal peak. The 
dependence of the width of the peak on different kinematic variables of the B^ decay 
was investigated. The strongest correlation was found to be with the angle between the 
momenta of the D*~ candidate and the bachelor vr"*" in the lab frame, denoted as ^bach- 
Simulated pseudo-experiments were used to find an optimal number of 6'bach bins to be 
used in a simultaneous fit. The outcome is that five bins are used, with ranges 0-0.046, 
0.046-0.067, 0.067-0.092, 0.092-0.128 and 0.128-0.4 rad, chosen to have approximately 
equal numbers of B^ decays in each. The peak width in the highest bin is approximately 
60 % of that in the lowest bin. The pseudo-experiments show that the simultaneous fit 
in bins of ^bach is approximately 20 % more sensitive to a potential B^ signal than the fit 
without binning. 

The signal yields are obtained from a maximum likelihood fit to the D*~7r^ invariant 
mass distribution in the range 5150-5600 MeV/c^. The fit is performed simultaneously in 
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the five 6'bach bins. Tlie fit includes double Gaussian shapes, where the two Gaussian func- 
tions share a common mean, for and signals, together with an exponential compo- 
nent for the partially reconstructed background, a linear component for the combinatorial 
background and a non-parametric function, derived from simulation, for B^ D*~K^ 
decays. The probability density function (PDF) for the B^ — )■ D*~K^ background is 
shifted by the mass difference between data and simulation for each bin of 6'bach- 

The parameters of the double Gaussian shapes are constrained to be identical for B^ 
and -Bg signals, with an offset in their mean values fixed to the known B^-B^ mass differ- 
ence [8] . Additionally, the relative normalisation of the two Gaussian functions and the 
ratio of their widths are constrained within uncertainties to the value obtained in simula- 
tion. A total of thirty-three parameters are allowed to vary in the fit: the ratio of yields 
N (B^) / N (B*^) , the linear slope of the combinatorial background and the exponential pa- 
rameter of the partially reconstructed background, plus separate parameters in each of 
the ^bach bins to describe the peak position and core Gaussian width of the signal PDF, 
and the yields of the B^ peak, the combinatorial background, the partially reconstructed 
background, and the background from B^ — )■ D*~K^. 

The results of the fit are shown in Fig. |2J The total number of B^ D*^7r^ decays 
is found to be 30 000 ±400, and the ratio of yields is determined to be N{B^)/N{B^) = 
(1.4 ±3.5) X 10"'^, where the uncertainty is statistical only. The number of B^ — t- D*~K^ 
decays found is 1 200 ± 200, with a correlation of 7% to the ratio of signal yields. 

The ratio of yields is converted to a branching fraction following 

B(Bl = ^ X g| X I X B{B» ^ D-n*) . (1) 

where e{B^) and €{B^) are the efficiencies for the B^ and B^ decay modes respectively, 
while fd (fs) is the probability that a b quark produced in the acceptance results in a 
[Bg) meson. Their ratio has been determined to be fs/ fd = 0.256 ± 0.020 [22] . 

The total efficiencies are (0.165 ± 0.002) % and (0.162 ± 0.002) % for the 5° and 5° 
decay modes respectively, including contributions from detector acceptance, selection cri- 
teria, PID and trigger effects. The ratio is consistent with unity, as expected. The PID 
efficiency is measured using a control sample of D*" — ^ Z)°7r~, K^tx" decays to 

obtain background-subtracted efficiency tables for kaons and pions as functions of their 
p and pt The kinematic properties of the tracks in signal decays are obtained from 
simulation, allowing the PID efficiency for each event to be obtained from the tables. Note 
that this calibration sample is dominated by promptly produced D* mesons. The remain- 
ing contributions to the total efficiency are determined from simulation, and validated 
using data. 

Systematic uncertainties on B(i?^ — > Z}*^7r^) are assigned due to the following sources, 
given in units of 1 x 10~^, summarised in Table [T] Event selection efficiencies for both 
modes are found to be consistent in simulation to within 2 %, yielding a systematic un- 
certainty of 0.02. The fit model is varied by replacing the double Gaussian signal shapes 
with double Crystal Ball [23] functions (with both upper and lower tails), changing the 
linear combinatorial background shape to quadratic and including a possible contribution 
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Figure 2: Simultaneous fit to tlie full data sample in five bins of ^bach^ (a) 0-0.046, (b) 0.046- 
0.067, (c) 0.067-0.092, (d) 0.092-0.128 and (e) 0.128-0.4 rad. Note the y-axis scale is logarithmic 
and is the same for each bin. Data points are shown in black, the full PDF as a solid blue line 
and the component PDFs as: (red dot-dashed) partially reconstructed background, (magenta 
dashed) combinatorial background, (blue dashed) signal, (black dot-dashed) signal and 
(green 3 dot-dashed) B^ — >■ D*~K~^ background. 
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Table 1: Systematic uncertainties on B{Bg D*^7r^). 



Source 


Uncertainty (10 ^) 


Efficiency 


0.02 


Fit model 


1.44 


Fit bias 


0.12 


Multiple candidates 


0.22 


fs/ fd 


0.12 




0.08 


Total 


1.47 



from D*~K^. The non-parametric function for the B'^ — > D*~K^ background 

was scaled in each bin to account for the change in the width of the B^ signal. Com- 
bined in quadrature these sources contribute 1.44 to the systematic uncertainty. Possible 
biases in the determination of the fit parameters are investigated by simulated pseudo- 
experiments, leading to an uncertainty of 0.12. Events with multiple candidates are 
investigated by performing a fit having chosen one candidate at random. This fit is per- 
formed 100 times, with different seeds, and the spread of the results, 0.22, is taken as 
the systematic uncertainty. The uncertainty on the quantity fs/fd contributes 0.12, while 
that on B{B'^ — > D*^7i^) gives 0.08. Combining all sources in quadrature, the total abso- 
lute systematic uncertainty is 1.47 x 10~^, and the B^ branching fraction is determined 
to be -> D*^n'^) = (1.5 ± 3.8 ± 1.5) x 10"^, where the first uncertainty is statistical 

and the second is systematic. 

A number of cross-checks are performed to test the stability of the result. Candidates 
are divided based upon the hardware trigger decision into three groups; events in which a 
particle from the signal decay created a large enough cluster in the calorimeter to fire the 
trigger, events that were triggered independently of the signal decay and those events that 
were triggered by both the signal decay and the rest of the event. The neural network 
and PID requirements are tightened and loosened. The non-parametric PDF used to 
describe the background from D*K decays was smoothed to eliminate potential statistical 
fiuctuations. All cross-checks give consistent results. 

Since no significant signal is observed, upper limits are set, at both 90 % and 95 % 
confidence level (CL), using a Bayesian approach. The statistical likelihood curve from 
the fit is convolved with a Gaussian function of width given by the systematic uncertainty, 
and the upper hmits are taken as the values containing 90 % (95 %) of the integral of the 
likelihood in the physical region. The obtained limits are 

-8(5° ^ D*^TT^) < 6.1 (7.8) X 10"^ at 90% (95%) CL . 

In summary, the decay i?^ D*^tx^ is searched for in a data sample of 1.0 fb~^ of 
data collected with the LHCb detector during 2011. No significant signal is observed 
and upper limits on the branching fraction are set. The absence of a detectable signal 
indicates that rescattering effects may make significant contributions to other hadronic 
decays, such as 5° vt+tt" and 5° DD, as recently suggested |5]. 
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